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Abstract 

A search for first-generation leptoquarks was performed in electron-proton and positron- 
proton collisions recorded with the ZEUS detector at HERA in 2003-2007 using an 
integrated luminosity of 366 pb _1 . Final states with an electron and jets or with miss- 
ing transverse momentum and jets were analysed, searching for resonances or other 
deviations from the Standard Model predictions. No evidence for any leptoquark 
signal was found. The data were combined with data previously taken at HERA, 
resulting in a total integrated luminosity of 498 pb" 1 . Limits on the Yukawa coup- 
ling, A, of leptoquarks were set as a function of the leptoquark mass for different 
leptoquark types within the Buchmuller-Riickl-Wyler model. Leptoquarks with a 
coupling A = 0.3 are excluded for masses up to 699 GeV. 
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1 Introduction 



Many extensions of the Standard Model (SM) predict the existence of particles carrying 
both baryon and lepton number, such as leptoquarks (LQs) pQ. In ep collisions at HERA, 
such states could have been produced directly through electroi]0-quark fusion (Fig. [T^i) if 
their masses, Mlq, were lower than the HERA centre-of-mass energy, yfs. The leptoquarks 
would have decayed into an electron and quark or an electron neutrino and quark, yielding 
peaks in the spectra of the final-state lepton-jetg^l invariant mass, My s . Leptoquarks with 
Mlq > y/s could not have been produced as resonances, but they would still have caused 
deviations from the SM prediction in the observed My s spectrum due to virtual LQ ex- 
change (Fig. [T^i and b). This paper presents an analysis of the My s spectrum searching 
for evidence for leptoquarks. 

The prediction for the M\j s spectrum is given by SM neutral current (NC) and charged 
current (CC) deep inelastic scattering (DIS) (Fig. [lb). Any leptoquark signal would have 
to be identified as a deviation from this SM prediction. At high Mij 8 , the SM prediction 
falls rapidly due to the dependence of the DIS cross sections on Q 2 , the virtuality of 
the exchanged boson, and to the sharply falling valence-quark density at large Bjorken 
x. This makes the data especially sensitive to virtual leptoquark exchange and LQ-DIS 
interference. 

Leptoquarks have been searched for previously in ep collisions [2H5] and in e + e~ [HIE], 
PP 0E] and pp pTJHT2] collisions. Using ep collisions, the Yukawa coupling, A, of pos- 
sible LQ states to electron and electron neutrino is probed. In pp and pp collisions, the 
LQ production proceeds via the strong interaction and is independent of A. Thus the 
experimental approaches are complementary and ep collisions provide a unique testing 
ground. For this paper, the predictions for LQ cross sections were determined in leading 
order (LO) using the CTEQ5D parton density functions |13| (PDFs) using the Buchmuller- 
Riickl-Wyler model [T|. This model assumes that some of the leptoquarks are doublets or 
triplets with degenerate masses. This differs from the assumptions made for production via 
the strong interaction where in general a singlet state is assumed. The LHC experiments 
so far provided only limits for scalar LQs [T0HT2] . 

In the analysis presented here, no evidence for any leptoquark signal was found. Therefore 
limits on A were derived as a function of M L q for the different leptoquark states described 
by the Buchmuller-Riickl-Wyler model. 

The analysis is based on the data collected by the ZEUS experiment in the period 2003- 
2007, corresponding to an integrated luminosity of 366 pb _1 . During this period, HERA 

1 Unless otherwise specified, 'electron' refers to both positron and electron and 'neutrino' refers to both 

neutrino and antineutrino. 
2 There can be more than one jet in the final state due to QCD initial or final state radiation. 
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was operated with a polarised lepton beam. The four data subsamples with different 
polarisation and lepton beam charge are summarised in Table [TJ The experimental set-up 
described in Section [2] is that corresponding to these subsamples. The final limits set also 
included data collected in 1994-2000, giving a total integrated luminosity of 498 pb _1 . 
Thus all data from ZEUS were included and the results presented here supersede those 
published previously [2H3]. 

2 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere |I4j. A brief outline 
of the components that are most relevant for this analysis is given below. 

In the kinematic range of the analysis, charged particles were tracked in the central tracking 
detector (CTD) |15| and the microvertex detector (MVD) j 16) . These components operated 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consisted of 72 cylindrical drift-chamber layers, organised in nine superlayers covering the 
polar-ang region 15° < 6 < 164°. The MVD silicon tracker consisted of a barrel (BMVD) 
and a forward (FMVD) section. The BMVD contained three layers and provided polar- 
angle coverage for tracks from 30° to 150°. The four-layer FMVD extended the polar-angle 
coverage in the forward region to 7°. After alignment, the single-hit resolution of the 
MVD was 24 ]im. The transverse distance of closest approach (DCA) to the nominal 
vertex in X—Y was measured to have a resolution, averaged over the azimuthal angle, of 
(46 © 122/pr) |^m, with p T in GeV. For CTD-MVD tracks that pass through all nine CTD 
superlayers, the momentum resolution was cr(p T )/p T = 0.0029j9y © 0.0081 © 0.0012/pr, 
with pt i n GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [17] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
was subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter was called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, were a(E)/E = 0.18/ \f~E for electrons 
and cr(E)/E = 0.35 /VE for hadrons, with E in GeV. 

The luminosity was measured using the Bethe-Heitler reaction ep — > eyp by a luminosity 
detector which consisted of independent lead-scintillator calorimeter [18] and magnetic 
spectrometer [19] systems. 

3 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
nominal proton beam direction, referred to as the "forward direction", and the X axis pointing towards 
the centre of HERA. The coordinate origin is at the centre of the CTD. The pseudorapidity is defined 
as rj = — In (tan |), where the polar angle, 9, is measured with respect to the Z axis. 
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The lepton beam in HERA became naturally transversely polarised through the Sokolov- 
Ternov effect [2"0"|[2T|. The characteristic build-up time for the HERA accelerator was 
approximately 40 minutes. Spin rotators on either side of the ZEUS detector changed 
the transverse polarisation of the beam into longitudinal polarisation and back again. The 
lepton beam polarisation was measured using two independent polarimeters, the transverse 
polarimeter (TPOL) [22] and the longitudinal polarimeter (LPOL) [23]. Both devices 
exploited the spin-dependent cross section for Compton scattering of circularly polarised 
photons off electrons to measure the beam polarisation. The luminosity and polarisation 
measurements were made over time intervals that were much shorter than the polarisation 
build-up time. The polarisation values measured in different data taking periods are shown 
in Table [TJ The fractional systematic uncertainty on the measured polarisation was 4.2% 
for TPOL and 3.6% for LPOL. 



3 Monte Carlo simulation 



Monte Carlo techniques were used to determine the SM DIS background and the Mij S 
resolution of a possible signal. 

Standard Model NC and CC DIS events were simulated using the Heracles 4.6.6 [21] 
program with the DjANGOH 1.6 [25] interfaces to the hadronisation programs and using 
CTEQ5D [13| PDFs. Radiative corrections for initial- and final-state electroweak radiation, 
vertex and propagator corrections, and two-boson exchange were included. The colour- 
dipole model of ARIADNE 4.12 [26J was used to simulate 0(ag) plus leading-logarithmic 
corrections to the result of the quark-parton model. ARIADNE uses the Lund string model 
of JETSET 7.4.1 |27| for the hadronisation. 

The production and decay of resonances were simulated using PYTHIA 6.1 [28] . which 
takes into account the finite width of the resonant state, but includes only the s-channel 
diagrams. It also takes into account initial- and final-state QCD radiation from the quark 
and the effect of LQ hadronisation before decay [29| as well as the initial-state QED 
radiation from the electron. Such simulated samples of LQ events were used to study the 
bias and the resolution for the reconstructed LQ mass. The prediction for a LQ contribution 
to the DIS samples was evaluated by reweighting the DIS samples according to the LQ 
production processes (Section [6]). 

The ZEUS detector response was simulated using a program based on GEANT 3.21 |30| . 
The generated events were passed through the detector simulation, subjected to the same 
trigger requirements as the data and processed by the same reconstruction programs. 
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4 Leptoquark signal expectation 



The Buchmuller-Riickl-Wyler (BRW) model [TJ was used to calculate the cross sections in 
LO using the CTEQ5D PDFs for the LQ states listed in Table [2j The next-to-leading-order 
(NLO) QCD corrections, the so-called i^-factors |31[l32]. available for scalar LQs were not 
considered, because equivalent calculations are not available for vector LQs. All limits 
presented in this paper are for LQ production in this model. 

As the final states for LQ production are identical to states produced in DIS, the DIS 
Monte Carlo samples were reweighted accordingly to produce predictions for the respect- 
ive Mij S spectra. The BRW model predicts a dependence of the cross sections on the 
beam polarisation. Therefore predictions were computed taking into account the average 
polarisation of the respective data samples. 

The BRW model includes both u and s channel and interferences with DIS. For Mlq > y/s, 
the full LQ cross sections were used. For the virtual exchange, the cross section has 
a A 4 dependence. The interference terms provide a contribution with a A 2 dependence. 
For Mlq < ^/s, some simplifications were introduced. While the interference terms were 
calculated as for large Mlq, the w-channel contribution, expected to be small, was neglected 
and the narrow-width approximation (NWA) was used for the resonant s-channel LQ 
production. 

The Born-level cross section for the resonant (s-channel) LQ production in the NWA is: 

a NWA = (J+l)^A 2 g(xo,M 2 Q )(l±P e ), 

where q(xo, M 2 q) is the initial-state quark (or antiquark) parton-density function in the 
proton for Xq = M 2 q/s, J is the spin of the LQ and the term 1 ± P e accounts for the 
dependence on the beam polarisation. In e~p (e + p) scattering, the polarisation dependence 
is given by 1 + P e (1 — P e ) for LQs coupling to right-handed fermions and 1 — P e (1 + P e ) 
for LQs coupling to left-handed fermions. 

The expected width of a LQ state with a mass between 100 and 300 GeV and a Yukawa 
coupling A = 0.1 ranges from 0.01 to 0.2 GeV. This justifies the use of the NWA for the 
s-channel contribution. To simulate a LQ signal in the MC, the s-channel term was added 
to the interference terms and the DIS MC events were reweighted. To reduce statistical 
fluctuations, the prediction from the NWA was smeared with a Gaussian narrower than 
the experimental resolution before adding the interference terms. 

The effect of QED initial-state radiation, which decreases the production cross section, was 
taken into account for both resonant and non-resonant LQ production. The effect is larger 
for resonant LQ production and ranges up to 25% for Mlq close to -y/s. 
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The polarisation dependence is expected to be different for LQ production and DIS. The 
separation of the data according to polarisation therefore provides a handle to identify a 
possible LQ signal. 

Another possibility to isolate a leptoquark signal is the angular dependence of the process. 
The variable 9*, the lepton scattering angle in the lepton-jets centre-of-mass frame, can be 
used to improve the signal-to-background ratio, especially for resonance production. The 
decay of a scalar resonance, for example, will result in a flat distribution in cos#*, while 
NC DIS events show approximately a 1/(1 — cos6**) 2 distribution. 



5 Signal search 

Events from a hypothetical resonance decaying into eq (uq) have a topology identical to 
DIS NC (CC) events. Hence the final state from a high-mass resonance is expected to 
have at least one jet and either an identified final-state electron or large missing transverse 
momentum. The lepton-jet invariant mass was calculated as 

M ljs = yj (1) 

where Eij S is the sum of the energies of the outgoing lepton and the selected jets and pij S is 
the vector sum of the lepton and jets momenta. The modulus of the transverse momentum, 
Pt, and the net transverse energy, E?, are defined as 



Pt = Px+Py= Ei sin 0i cos & j + (^Yl 
E T = 22 Ei sin 9i, 



2 

Ei sin 6i sin fa I , (2) 



where the sum runs over all calorimeter energy deposits, E t . The polar and azimuthal 
angles, 6^ and fa, of the calorimeter energy deposits were measured relative to the recon- 
structed event vertex. The quantity E — pz, also used in the event selection, is defined 

as 

E - p z = ^2(E-p z ) i = ^2(E i -E i cob 9i). (3) 

i i 

The hadronic polar angle was calculated as |33| 

^T.had + {E- ^Zjhad 
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where -P^had and {P — Pz) had are calculated as in Equations (j2j) and ([3]), but with the sum 
running only over the calorimeter energy deposits belonging to the hadronic final state. In 
case of the CC topology, P^taA an d {P ~ -Pz)had are equivalent to p\ and E — pz from 
Equations (j2J) and (j3J), respectively. 

5.1 Neutral current, ep — >- eX, topology 
Event selection 

The double-angle (DA) method [M] was used to reconstruct the kinematic variables. In 
this method, the polar angle of the scattered electron and the hadronic polar angle (see Eq. 
PJ) were used to reconstruct the kinematic variables £da ; Vba, and Q^A- The inelasticity 
y was also reconstructed as y e , using the electron method |34] . 

Events with the topology ep — » eX, where X denotes one or more jets, were selected using 
the following criteria: 

• the Z coordinate of the reconstructed event vertex was required to be in the range 
\Z\ < 30 cm, consistent with an ep collision; 

• an electron identified in the CAL [35] was required. If the electron was found within 
the acceptance of the tracking detectors, a track matched to the energy deposit in 
the calorimeter was required with the distance of closest approach between the track 
extrapolated to the calorimeter surface and the energy cluster position to be less 
than 10 cm and the electron track momentum, p^ k , to be larger than 3 GeV. A 
matched track was not required if the electron emerged at a polar angle outside the 
acceptance of the tracking detector. Instead the electron was required to have a 
transverse momentum greater than 30 GeV. An isolation requirement was imposed 
such that the energy not associated with the electron in an r\ — <fi cone of radius 0.8 
centred on the electron was less than 5 GeV; 

• a fiducial- volume cut was applied to the electron to guarantee that the experimental 
acceptance was well understood. It excluded the transition regions between the FCAL 
and the BCAL. It also excluded the regions within 1.5 cm of the module gaps in the 
BCAL; 

• at least one hadronic jet with transverse momentum p 3 T > 15 GeV and \rj\ < 3, 
obtained using the kx cluster algorithm [36] in its longitudinally invariant inclusive 
mode [37], was required. The centroid of any jet at the FCAL face was required 
to be outside a box of 40 x 40 cm 2 centred on the proton beam [38], in order to 
ensure good energy containment and to reduce the systematic uncertainties due to 
the proton remnant. Additional jets were required to have p 3 T > 10 GeV and \r]\ < 3; 
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• to restrict the phase space to the region most relevant to the LQ search, cuts on 
Ql A > 2500 GeV 2 and x DA > 0.1 were applied. Higher values of x and Q 2 correspond 
to higher LQ masses, where the signal-to-background ratio for leptoquark events is 
higher; 

• to avoid phase-space regions in which the DIS MC generator was not valid, the 
quantity ?/da(1 — ^da) 2 was required to be larger than 0.004; 

• to remove background from phot opro duct ion events and beam-gas events overlaid on 
NC events, the requirements 38 < E — pz < 65 GeV and y e < 0.95 were imposed. To 
remove cosmic-ray events and beam- related background events, the quantity Pt/V E? 
was required to be less than 4V GeV and the quantity p?/ E? was required to be less 
than 0.7. Elastic QED Compton events were rejected by considering the balance 
between the photon and the electron [39 1. 

The mass shifts and resolutions for resonant lepton-quark states were calculated from 
the LQ MC. The mass resolution, determined from a Gaussian fit to the peak of the 
reconstructed mass spectrum, fell from 5% to 3% as the resonant mass increased from 150 
to 290 GeV. Any mass shift was within 0.5% for LQ masses between 150 and 290 GeV. 



Search results 

After the above selection, 9 369 events were found in the data from 2003-2007, compared to 
9 465 ± 494 expected from the NC MC and the evaluation of its systematic uncertainties 
(see below). The measured distributions of the M e j s spectra for e~p (e + p) data with a 
left-handed and a right-handed lepton beam are shown in Figs. [2] and [3] (Figs. [4] and [5]), 
respectively. A cut on cos 8* < 0.4, introduced to suppress further the SM background [2jl3] . 
was also applied. The upper parts of the plots show the spectra with and without the cos 8* 
cut, while the lower parts show the ratio of the observed spectrum to SM expectations 
with no cut applied on cos 8*. Good agreement is seen between the data and the SM 
NC prediction. Figures [6] and [7] show the e~p data for the left-handed electron beam 
together with the predictions for a Sq LQ state with a mass of 210 GeV and a coupling A 
of 0.3 as well as a mass of 400 GeV and a coupling A of 1 without and with a cos 8* cut, 
respectively. 



Systematic uncertainties 

The uncertainty on the expected number of events from SM NC DIS processes was invest- 
igated. Relevant were [33J: 
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• the uncertainty on the calorimeter energy scale, 1% for electrons and 2% for hadrons. 
This led to an uncertainty of 1% (6%) in the NC expectation for M e j S =150 (220) GeV; 

• the uncertainty on the parton densities as estimated by Botje [ID], which gave an 
uncertainty of 3% (5%) in the NC expectation for M ejs =150 (220) GeV; 

• the uncertainty on the luminosity determination of 1.8% for electron and 2.2% for 
positron beamd3, which is directly reflected in the result. 

The overall systematic uncertainties on the background expectations were obtained by 
adding all relevant contributions in quadrature. They are shown as hatched bands in 
Figs. [21 El SI and |5j For a given mass, the systematic uncertainty for the LQ signal is 
assumed to be the same as for the SM background. 

5.2 Charged current, ep — >• uX, topology 
Event selection 

The events with the topology ep — > vX, where X denotes one or more jets, are classified 
according to 7 , the hadronic polar angle (see Eq. (T4])) assuming a nominal vertex position 
of Z = 0. Events for which the hadronic system is not contained in the CTD acceptance 
(7o < 0.4 rad) are called low-70 events. The hadronic systems of high-70 events with 
7o > 0.4 rad are inside the CTD acceptance, so that cuts based on tracking information 
can be applied. The kinematic variables were reconstructed using the Jacquet-Blondel 
method [41] . 

The events were selected using the following criteria: 

• the Z coordinate of the reconstructed event vertex was required to be in the range 
\Z\ < 30 cm, consistent with an ep collision; 

• to restrict the phase space to the region most relevant to the LQ search, a cut on 
Qj B > 700 GeV 2 was applied. Since the resolution on Qj B was poor at high y, a cut 
on ?/jb < 0.9 was added; 

• a missing transverse momentum px > 22 GeV was required and, to suppress beam- 
gas events, the missing transverse momentum excluding the calorimeter cells adjacent 
to the forward beam hole was required to exceed 20 GeV. Compared to the ZEUS 
CC DIS analyses [421444] . these cuts are more stringent because CC events with low 
Pt lead to low invariant masses of the LQs and are therefore not of interest in this 
analysis; 

For a fraction of the positron data, the uncertainty was 3.5%, while for most of the positron data the 
uncertainty was 1.8%. This lead to a total uncertainty of 2.2%. 



8 



• in the high-70 region, cuts based on the number and quality of tracks were applied. 
Tracks with a transverse momentum above 0.2 GeV were selected. They were required 
to start from the MVD or the innermost superlayer of the CTD, and had to reach 
at least the third superlayer of the CTD. If in addition they pointed to the primary 
vertex, they were considered as "good tracks". At least one good track was required 
in the event and the ratio of the total number of tracks to the number of good tracks 
had to fulfill the requirement Nt * k ~ 20 < 4, where A^rk is the total number of tracks 

* ''gtrk 

and -/Vgtrk is the number of good tracks. This cut removed beam-gas events which are 
characterised by a high number of poor-quality tracks; 

• at least one hadronic jet with transverse momentum p 3 T > 10 GeV and \r]\ < 3, 
obtained using the cluster algorithm [36] in its longitudinally invariant inclusive 
mode [37], was required. The centroids of all jets at the FCAL face were required 
to be outside a box of 40 x 40 cm 2 centred on the proton beam [38], in order to 
ensure good energy containment and to reduce the systematic uncertainties due to 
the proton remnant. 

• to reject photoproduction and di-lepton background, for events with p? < 30 GeV a 
dedicated cut based on the energy distribution in the detector relative to the total 
transverse momentum was used. The transverse momentum sum for the calorimeter 
cells with a positive contribution to the total transverse momentum (parallel com- 
ponent Vp) and the corresponding sum for cells giving a negative contribution to 
the total transverse momentum (antiparallel component Vap) had to satisfy the con- 
dition Vap/Vp < 0.35 jll]. This requirement demanded an azimuthally collimated 
energy flow. In addition, for events with at least one good track, the azimuthal-angle 
difference, A0, between the missing transverse momentum measured by the tracks^) 
and that measured by the calorimeter was required to be less than 45° for events 
with px < 30 GeV and less than 60° otherwise. This cut rejects events caused by 
cosmic rays or muons in the beam halo as well as beam-gas events; 

• NC events were removed by discarding events containing electron candidates with an 
energy greater than 4 GeV [43j; 

• requirements on energy fractions in the calorimeter cells plus muon-finding algorithms 
based on tracking, calorimeter and muon-chamber information were used to reject 
events caused by cosmic rays or muons in the beam halo. Furthermore, the deposition 
times of the energy clusters in the calorimeter were checked to be consistent with 
the bunch- crossing time to reject events due to interactions between the beams and 
residual gas in the beam pipe or upstream accelerator components. In addition, 
topological cuts on the transverse and longitudinal shower shape were imposed to 
reject beam-halo muon events that produced a shower inside the FCAL. Cuts on 

5 The missing transverse momentum measured by the tracks is calculated using all the good tracks. 
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the calorimeter cell with the highest transverse energy were applied to reject sparks 
faking a CC event [33J. 

The neutrino energy and angle were calculated by assuming that missing and missing 
E — Pz were carried away by a single neutrino and used to calculate the invariant mass of 
the z/-jets system, M„j 8 , according to Eq. (JTJ). 

The shift and resolution of the invariant mass were studied by using the LQ MC events 
and fitting the mass peak with a Gaussian function. The resulting mass shift was within 
0.5% for LQ masses between 150 and 290 GeV, with the resolution varying from 8% to 6%, 
respectively. 



Search results 

After the above selection, 8 990 events were found in the data from 2003-2007, compared 
to 9 068 ± 501 expected from the CC MC and the evaluation of its systematic uncertainties 
(see below). The measured distributions of the M V j S spectra for the left-handed and right- 
handed e~p (e + p) data are shown in Figs. [S]and|H] (Figs. [TU1 and [TT]) . The upper parts of 
the plots show the spectra with and without the cut cos#* < 0.4, while the lower parts 
show the ratio of the observed spectrum to SM expectations with no cut applied on cos#*. 
Good agreement is seen between the data and the SM CC prediction. 



Systematic uncertainties 

The uncertainty on the predicted background from SM CC DIS processes was investigated. 
The uncertainties found to be relevant [33] are similar to those described in Section I5TT1 for 
the ep — > eX case and arise from: 

• the uncertainy on the hadronic energy scale of 2%, which led to an uncertainty of 
3% (10%) in the CC expectation for M vjs =150 (220) GeV; 

• the uncertainty on the parton densities as estimated by Botje (40] , giving 3% (4%) 
and 7% (9%) uncertainties on the cross section for e~p and e + p, respectively, for 
M U j S =150 (220) GeV. The correlations between e~p and e + p as well as NC and CC 
cross section uncertainties were taken into account; 

• the uncertainty on the luminosity determination of 1.8% for electron and 2.2% for 
positron beams, which is directly reflected in the result; 
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• the uncertainty on the measured polarisation of 4.2%. To be conservative, the TPOL 
uncertainty, which is larger than the LPOL uncertainty, was used. This led to an 
uncertainty on the SM cross section of 0.9% (2.4%) for left-handed e~p (e + p) data 
and 1.8% (1.0%) for right-handed e~p (e + p) data, respectively. 

The overall systematic uncertainties on the background expectations were obtained by 
adding all relevant contributions in quadrature. They are shown as hatched bands in 
Figs. [HI El EH] and [TTJ For a given mass, the systematic uncertainty for the LQ signal is 
assumed to be the same as for the SM background. 



6 Limits on leptoquarks 



The expectation from a potential LQ signal was obtained by reweighting the DIS MC 
according to the cross sections predicted in the BRW model (see Section H]). Each MC 
event is reweighted with the following weighting factor, WF: 



WF(x, y- P e - M LQ , A) = dxdy J M ' 



dxdy 



where x and y are the true kinematic variables of the MC simulation, and P e is the average 
polarisation of the data sample given in Table [TJ The effect of QED initial-state radiation 
was taken into account. 

The limits were calculated including the results of the search presented here and the data 
recorded with the ZEUS detector in the years 1994-2000 jl]. They were set using a binned 
likelihood technique in the (My s , cos 6*) plane. The region 150 < M; JS < 320 GeV was 
used. The data were binned separately for each of the data sets listed in Table [TJ thereby 
taking into account different beam charges and polarisation. For leptoquark states with 
vq decays, both the eq — > eX and the eq — > vX samples were used, while for leptoquark 
states decaying only to eg, only the eq — > eX samples were used. 

The upper limit on the coupling strength, Aumit, as a function of Mlq, was obtained by 
solving 

A limit OO 

J dX 2 L(M LQ ,A) = 0.95 J d\ 2 L(M LQ ,X), (5) 
o o 

where L is the product of the Poisson probabilities of the observed number of events in 
bin i, Ni, with the expectation, fii, convoluted with Gaussian distributions for the main 
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systematic uncertainties, yielding a modified expectation fi'f. 



oo 

-oo i v i 

The index j denotes the source of a systematic uncertainty and 5j corresponds to the 
variation of the j th systematic parameter in units of the nominal values quoted in Sections 
15.11 and 15.21 The index i runs over the bins in My s -cos 6* and the relevant data sets. The 
modified expectation was calculated as 

j 

where gives the fractional variation of /Xj under the nominal shift in the j th systematic 
parameter. This ansatz of ii- reduces to a linear dependence of fi' { on each 5j when 6j is 
small, while avoiding the possibility of /i- becoming negative which would arise if ^ was 
defined linear function of the 8j. 

The coupling limits for the 14 BRW LQs listed in Table [2] were calculated for masses up to 
1 TeV. For large LQ masses, i.e. Mlq 3> y/s, their exchange can be described as a contact 
interaction with an effective coupling proportional to X 2 /MIq. Table [3] shows the limits on 
A for all BRW LQs at M LQ = 1 TeV. 

Figures [T2Hl~5l show the coupling limits on the scalar and vector LQs with F = and 
F = 2, respectively, where F = 3B + L is the fermion number of the LQ and B and L 
are the baryon and lepton numbers, respectively. The limits range from 0.004-0.017 for 
Mlq = 150 GeV, and from 0.43-3.24 for Mlq = 1 TeV. The lowest masses for which LQs 
with A = 0.1 and with A = 0.3 are not excluded are summarised in Table HI They range 
from 274 to 300 GeV for A = 0.1 and from 290 to 699 GeV for A = 0.3. 

The limit on the LQ state Sy 2 (Sq) can be interpreted as a limit on \\f]3 for an up-type 

squark ul (a down-type squark dp) in supersymmetric models with i?-parity violation |45) . 
where A is the coupling of ul to eq (dp to eq and vq) and /3 is the branching fraction of 
the squarks to lepton (e or v) and quark£l 

Figure [16] f lTTl) shows the limits on the S^ 2 (Sf) LQ compared to the limits from AT- 
LAS HO), HI |5| L3 \6\ and OPAL [7J. The L3 and OPAL limits were evaluated up to 
500 GeV. Limits using pp or pp collisions are obtained from leptoquark pair production, 
which is independent of A. Therefore it is not obvious whether the limits should be com- 
pared directly. The ATLAS exclusion range given in Fig. [17] also depends on the assump- 
tion that the branching ratio of the LQ state to electron and quark is one. Limits using 

6 The branching fractions of the squarks to eq and vq are assumed to be j3 eq = ft, j3 uq = for ul, and 
I3 eq = 0.5/3, p uq = 0.5/3 for d R . 
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e + e~ collisions are obtained from a search for indirect effects in the process e + e~ — > qq. 
In general, the limits from this analysis are significantly better than the LEP limits for 
Mlq < y/s, and comparable for Mlq > \/s. The limits obtained by ZEUS are similar to 
those obtained by HI. 



7 Conclusions 

Data recorded by the ZEUS experiment at HERA were used to search for the presence 
of first-generation scalar and vector leptoquarks. The data samples include 185 pb" 1 of 
e~p and 181 pb~ of e + p collisions with polarised electrons and positrons. No resonances 
or other deviations from the Standard Model were found. The inclusion of data with 
unpolarised beams yields a total set of data corresponding to 498 pb" 1 , which was used to 
set upper limits on the Yukawa coupling A for the 14 Buchmuller-Ruckl-Wyler leptoquark 
states as a function of the leptoquark mass. Assuming A = 0.3, the mass limits range from 
290 to 699 GeV. 
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period 


lepton 


luminosity (pb 1 ) 




(GeV) 


04-06 


e~ 


106 


-0.27 


318 


04-06 


e~ 


79 


0.30 


318 


03-04 




17 






06-07 


e+ 


60 


-0.37 


318 




03-04 




21 






06-07 


e+ 




0.32 


318 




83 


94-97 


e+ 


49 





300 


98-99 


e~ 


17 





318 


99-00 


e+ 


66 





318 



Table 1: Details, including longitudinal polarisation, P e , of the different data samples 
used. 
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LQ species 


charge 


production 


decay 


branching ratio 


coupling 
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1 




qR 
°l/2 


5/3 


etlir 


e + u 


1 


Ap 




2/3 


etdi 


e + d 


1 


-A fl 


CL 
°l/2 


-2/3 


e R d R 


e + d 


1 


Al 


V L 


2/3 
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Al 
Al 
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1 


\r 
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V? 
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e R U L 
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1 
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2/3 
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1/2 


-A L 
Al 


F=2 




oL 
D 


-1/3 


e L u L 


e u 
v e d 


1/2 
1/2 


Al 
-Al 


qR 
°0 


-1/3 


e R u R 


e~u 


1 


\r 


qR 
D o 


-4/3 


e~ R d R 


e~d 


1 


Xr 


si 


-1/3 


e~ L u L 


e~u 
v e d 


1/2 
1/2 


-Al 
-A l 




-4/3 


eld L 


e~d 


1 


-V2\ L 


Vl/2 


-4/3 


z~id R 


e~d 


1 


Al 


yR 


-4/3 


e~ R d L 


e~d 


1 


Xr 




-1/3 


e R u L 


e~u 


1 


x R 


V L 

V l/2 


-1/3 


elu R 


e~u 


1 


Al 



Table 2: Leptoquark species for fermion number F = and F = 2 as defined in the 
Buchmuller-Ruckl-Wyler model and the corresponding couplings. Those LQs that couple 
only to neutrinos and quarks and therefore could not be produced at HERA are not listed. 
The LQ species are classified according to their spin (S for scalar and V for vector), their 
chirality (L or R) and their weak isospin (0,1/2,1,). The leptoquarks S and V differ by 
two units of hypercharge from S and V , respectively. In addition, the electric charge of 
the leptoquarks, the production channel, as well as their allowed decay channels assuming 
lepton- flavour conservation are displayed. The nomenclature follows the Aachen convention 
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LQ type (F=0) 


V L 


v * 


yR 




oL 
°l/2 


cR 
°l/2 


°l/2 


Aiimit 


0.87 


1.91 


0.76 


0.43 


1.00 


2.29 


1.91 


LQ type (F=2) 


nL 


qR 
°0 


cR 
°0 


st 


*l/2 


V l/2 


V L 


A limit 


1.15 


1.48 


3.24 


0.60 


1.95 


0.95 


0.76 



Table 3: Upper limit on the Yukawa coupling, Aii m ; t as defined in Eq. |3]), for the 14 BRW 
LQs at M LQ = 1 TeV. 



LQ type (F=0) 




v R 




V{ 


qL 
°l/2 


ci?, 

°l/2 


oL 
°l/2 


M LQ (GeV) (A limit = 0.1) 


276 


275 


295 


300 


295 


294 


274 


M LQ (GeV) (AiMt = 0.3) 


325 


292 


376 


699 


345 


300 


292 


LQ type (F=2) 


oL 
°0 


qR 


qR 




V L 


K l/2 


V L 


M LQ (GeV) (A limit = 0.1) 


295 


292 


274 


298 


278 


293 


293 


M LQ (GeV) (A limit = 0.3) 


332 


304 


290 


506 


292 


303 


376 



Table 4-' Lower limit for the masses of the 14 BRW LQs for Aii m ; t =0.1 and \u m it=0.3 as 
deduced from Eq. |3p. 



20 




Figure 1: Diagrams for (a) s-channel LQ production/ exchange and (b) u-channel LQ 
exchange and for (c) SM deep inelastic scattering via photon, Z° and W exchange. 
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Figure 2: Comparison of the left-handed e~p sample (dots) and the NC SM expectation 
(solid histogram) for the reconstructed invariant mass, M e j S; in the e~p — > e~X topology. 
The data (open squares) and the SM expectation (dashed histogram) for cos 6* <0.4 are 
also shown. The shaded area shows the overall uncertainty of the SM MC expectation. The 
lower part of the plot shows the ratio between the data and the SM expectation without the 
cos 6* cut. 
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Figure 3: Comparison of the right-handed e~p sample (dots) and the NC SM expectation 
(solid histogram) for the reconstructed invariant mass, M e j S; in the e~p — > e~X topology. 
Other details as in the caption to Fig. 
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Figure 4 : Comparison of the left-handed e + p sample (dots) and the NC SM expectation 
(solid histogram) for the reconstructed invariant mass, M e j S , in the e + p — > e + X topology. 
Other details as in the caption to Fig. 
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Figure 5: Comparison of the right-handed e + p sample ( dots ) and the NC SM expectation 
(solid histogram) for the reconstructed invariant mass, M e j s , in the e + p — > e + X topology. 
Other details as in the caption to Fig. [H 
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Figure 6: Comparison of the reconstructed invariant mass, M e j S; distribution in the 
e~p — > e~X topology for the left-handed e~p sample (dots) to the NC SM expectation 
(solid histogram) and to the predictions of the model including a Sq LQ state with a mass 
of 210 GeV and a coupling A of 0.3 (dashed histogram) as well as a mass of 400 GeV and 
a coupling X of 1 (dotted histogram) . 
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Figure 7: As Fig. 0, but the cut on cos 6* <0.4 was applied. 
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Figure 8: Comparison of the left-handed e~p sample (dots) and the CC SM expectation 
(solid histogram) for the reconstructed invariant mass, M u j S; in the e~p — > vX topology. 
The data (open squares) and the SM expectation (dashed histogram) for cos 6* <0.4 are 
also shown. The shaded area shows the overall uncertainty of the SM MC expectation. The 
lower part of the plot shows the ratio between the data and the SM expectation without the 
cos 6* cut. 
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Figure 9: Comparison of the right-handed e~p sample (dots) and the CC SM expectation 
(solid histogram) for the reconstructed invariant mass, M U j S; in the e~p — > vX topology. 
Other details as in the caption to Fig. O 
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Figure 10: Comparison of the left-handed e + p sample (dots) and the CC SM expectation 
(solid histogram) for the reconstructed invariant mass, M V j S , in the e + p — > vX topology. 
Other details as in the caption to Fig. 0. 
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Figure 11: Comparison of the right-handed e + p sample (dots) and the CC SM expectation 
(solid histogram) for the reconstructed invariant mass, M U j S , in the e + p — > vX topology. 
Other details as in the caption to Fig. 0. 
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Figure 12: Coupling limits, A iim it, as a function of LQ mass for scalar F=0 BRW LQs. 
The areas above the curves are excluded according to Eq. (GJ). 
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Figure 13: Coupling limits, Aii m ; t , as a function of LQ mass for vector F=0 BRW LQs. 
The areas above the curves are excluded according to Eq. |3J). 



27 



ZEUS 



- 1 1 1 1 

F=2 scalar LQ limit 


i i i i ! i i i 




-ZEUS e*p (498 pb" 1 ) 






: & .-"^ 

& f 

: # 

- % L 




— ~ 


& ' i 

— ' i 

- 1 ' f 


— S L 




_ it 
it 
•* * / 


--S R 










mIimiIiimImiiImm 


st 

1 , , , , 1 , , . , 1 , , . , 1 





0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



M LQ (TeV) 



Figure 14: Coupling limits, Aiimit, as a function of LQ mass for scalar F=2 BRW LQs. 
The areas above the curves are excluded according to Eq. (GJ). 




Figure 15: Coupling limits, \ \\ m \t, as a function of LQ mass for vector F=2 BRW LQs. 
The areas above the curves are excluded according to Eq. |3]). 
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Figure 16: Coupling limits as a function of LQ mass for the Sy 2 LQ from ATLAS, L3, 
HI and ZEUS. 




Figure 17: Coupling limits as a function of LQ mass for the S± LQ from ATLAS, OPAL, 
HI and ZEUS. 
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